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ABSTRACT 



(N 

c 2 ^ ' Context. Eruption of a coronal mass ejection (CME) drags and "opens" the coronal 



magnetic field, presumably leading to the formation of a large-scale current sheet and 
the field relaxation by magnetic reconnection. 

Aims. We analyze physical characteristics of ray-like coronal features formed in the 
aftermath of CMEs, to check if the interpretation of this phenomenon in terms of 
reconnecting current sheet is consistent with the observations. 

Methods. The study is focused on measurements of the ray width, density excess, and 
coronal velocity field as a function of the radial distance. 

Results. The morphology of rays indicates that they occur as a consequence of Petschek- 
like reconnection in the large scale current sheet formed in the wake of CME. The 
hypothesis is supported by the flow pattern, often showing outflows along the ray, 
and sometimes also inflows into the ray. The inferred inflow velocities range from 3 to 
30 km s~^ , consistent with the narrow opening-angle of rays, adding up to a few degrees. 
The density of rays is an order of magnitude larger than in the ambient corona. The 
density-excess measurements are compared with the results of the analytical model in 
which the Petschek-like reconnection geometry is applied to the vertical current sheet, 
taking into account the decrease of the external coronal density and magnetic fleld 
with height. 

Conclusions. The model results are consistent with the observations, revealing that the 
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main cause of the density excess in rays is a transport of the dense plasma from lower 
to larger heights by the reconnection outflow. 

Key words. Sun: coronal mass ejections (CMEs) - Sun: corona - (Sun:) solar wind - 
magnetohydrodynamics (MHD) 



1. Introduction 

According to the current comprehension of solar coronal mass ejections (CMEs), the erup- 
tion of an unstable magnetic structure is tightly associated with the formation of the large 
scale current sheet (hereinafter CS) in the wake of the eruption. This concept, connect- 
ing the large scale eruption and the localized ener gy release in a form of fl are, w as put 
forward and initia lly developed by Carmichael ll964l : Sturrock ll966t Hiravama ll974i Kopp 
& Pneuman 1 19761 (the so-called CSHKP model). When the current sheet becomes long 



enough, the tearing instability sets in, leading to fast reconnection of magnetic field (for 
the analytical, numerical, laboratory and observational re sults see, e.g.. 



Furth et al 



1963 



Ugai 



1987 



Gekelman fc Pfister 



1988 



Vrsnak et al 



20031 respectively). Reconnection in 



the post-CME current sheet results in an abrupt energy release causing a flare, and on the 



Lin 


2004; 


Vrsnak 


2008 



ture (e.g.. iLinI 12004 IVrsnakll2008l and references therein). Nimrerical MHD simulations 
covering various scales, from low corona up to 1 AU, also demonstrated post-CME cur 
rent sheet formation, clearly revealing the importa nce of the reconnection in the initiation 



Tor ok et al. 


2004; 


Rilev et al. 


2007 



Rilev et al 



2002 : 



Roussev et al 



2003 



20071 for a review see Forbes et al. 2006). 
The most prominent consequence of CME-associated reconnection is the appearance of 
the so called two-ri bbon flare, closely synchronized with the CME acceleration stage (e.g.. 



Maricic et al 



20071 and references therein). The ribbon expansion away from the magnetic 
inversion line, associated with growing system of hot loops connecting the ribbons, has led 
to the form ulation of the CSHKP model. The discovery of cusped stru cture above the flare 
loops (e.g., 



Tsuneta et al 



sourc es (e.g.. 



1992' 



Masuda et al 



Forbes fc Acton 



1994 



2004 : 



2006 ) , the loop shrinkag e (e.g.. 



1996 



Aschwanden et al 



Svestka et al 



1987 



Tsuneta 



1996; 



19961) , loop t op hard X-ray 



Sui et al 



Forbes fc Acton 



2004 



Vrsnak et al. 



1996 ; 



Veronig et al 



Sheelev et al 



20061 ) ■ and the recognition of growing posteruption loop systems in the 
absence of flares, further supported the CSHKP scenario. Additional evidence for recon- 
necti on below the erupting flux-rope may be found in the so-calle d disconnection events 
(e-g-, 



Webb fc Cliver 



1995 



Simnett et al 



1997 



Wang et al 



McKenzie fc Hudson 



1998; 



199£), downflows above the 



Innes et al 



post-eruption arcades (e.g., 
horizontal converging flows above the loops (jYokovama et al 



2003' 



2001 



flare-associated radio emission below the eruptive prominence ( Vrsnak et al 



Recently, post-CME features seen in UV spectra ( Ciaravella et al 



Innes et al 



2003b 



X-rav images (ISiii et al 



Webb et al 



Lin et al 



2003 ; 



2005 



2004 



Lin et al 



Bemporad et al 



Asai et al 



Lin et al 



20031) . 



2004J) 



20051), and 



2006 : 



2002 



Ko et al 



Ciaravella fc Raymond 2008 ) 



20051) . and white light coronagraph images (|Ko et al 



2003 



2003 



20051 ) have been attributed to the current sheets expected to 
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"connect" the flare loops to the CME core. The UV spectral signature is generally emis- 
sion in the high temperature lines of Fc XVIII or FeXIX formed at 5-10 MK. EUV images 
from the Extreme-ultraviolet Imaging Telescope (EIT) on board the Solar and Heliospheric 
Observatory (SoHO) and the Transition Regi on and Coronal Explorer (TRACE) show emis - 
sion in FeXXIV at even higher temperatures ( McKenzie fc Hudsonlll999 : Innes et alJboOS ). 



The UV features decay slowly, on time scales f rom several hou rs (jCiaravella fc Raymond 



20081) up to a few days (jBemporad et al 



2006 



Ko et al 



20031 ). The X-ray data from the 



Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) show regions of high 
temperature emission above loop tops, and both the morphology at different energies and 
the evolution of the structures supports the CS indentification. The white light features 
are bright, narrow rays th at appear to map from the cusped flaring loop system to the 



CME core ( Lin et al 



20051 ). Some show blobs of plasma moving outwards at several hun- 



dred kms^^. Like the UV features, the white light structures can last for a day or two. In 
general, white light observations reveal the electron column density, UV data show tem- 
peratures, emission measures and Doppler shifts, and X-ray observations reveal very high 
temperature plasma. 

The aim of this paper is to demonstrate that post-CME rays appear as a consequence 
of the reconnection in the current sheet formed in the wake of CME, and to quantify their 
basic physical characteristics. We outline the working hypothesis in Sect. 2, whereas details 
of the model used to compare the theoretical expectations with observations are provided 
in the Appendix. Observations and measurements are presented in Sect. 3, focusing on 
the morphology and plasma densities in post-CME rays. The temperature structure and 
the ionization state will be presented in a separate paper. Our density measurements are 
compared with the model estimates and some previous empirical results in Sect. 4. Results 
are discussed and conclusions drawn in Sect. 5. 



2. Working hypothesis 



Since post-CME rays extend outward from cusped structures associated with growing post- 
eruption loop-systems, our working hypothesis is that they appear as a consequence of 
the reconnection in the CS formed in the wake of CME. Furthermore, we assume the 
reconne ction process c an be approximately described in terms of the steady-state Petschek 
regime (jpetscheklllQel . In this regime, the reconnection takes place within a small diffusion 
region (DR), out of which two pairs of standing slow- mode shocks (SMSs) extend along 
of the axis of s ymmetry. Thus, the CS is bifurcated, w ith SMSs separating the inflow and 



outflow region ([Petschek 



1964 : 



Soward fc Priest 



19821) . At SMSs the inflowing plasma is 



compressed, heated, and accelerated, forming the upward and downward reconnection jet 
(Fig- HI)- Note that electric currents are concentrated only in DR and SMSs. 

When reconnection takes place in the homogeneous environment, the characteristics of 
outflowing plasma are determined primarily by the external plasma-to-magnetic pressure 
ratio, wherea s the outflow speed is approximately eq ual to the external Alfven speed (see 
Appendix in 



Aurass et al 



2002 : 



Skender et al 



20031 ). However, in the case of a vertical 



current sheet, the ambient coronal density and magnetic field decrease with height, so 
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Fig. 1. a) Interpretation of the ray in terms of the iPetschelJ (jl964 ) reconnection modeh 
DR - diffusion region; SMS - slow mode shocks (dashed hnes intersecting at DR); PEL 

- post-eruption loops. Magnetic field-lines are drawn by thin lines; gray area between the 
SMSs outlines the post-CME ray. The coordinate system is indicated (the line of sight is 
in y-direction). b) An element of bifurcated current sheet; arrows mark plasma flows (uin 

- inflow velocity; fgw ~ solar wind speed; Wout ~ outflow speed), where the arrow thickness 
depicts plasma densities. 
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Fig. 2. Dependence of the width and morphology of the ray at a given height on the ori- 
entation/geometry of the current sheet: a) straight CS; b) wavy CS. The CS half-thickness 
is denoted by d, horizontal length (i.e., the length in a;?/-plane, see also Fig. I) by L, the 
ray width by Z?, and the inclination to the line of sight (LoS) by 0. The plasma-column 
length A is drawn by thick double-arrow. At the bottom the intensity profile is sketched. 



the characteristics of the reconnection jets depend also on the height. The main effect is 
transport of the dense plasma from lower heights upward, making the outflow jet much 
denser than the ambient corona (see Appendix). In this respect, the post-CME ray should 
be, to a certain degree, similar to coronal streamers. 
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Bearing in mind the geometry of the Petschek-type reconnection, the post-CME CS 
should show distinct morphological characteristics. It should be thinnest at the height 
of the diffusion region, getting wider with increasing height. In the hom ogeneous plasma 
the angle between the SMSs is determined by the inflow Mach number (jSoward fc Priest 



1982: 



Vrsnak fc Skender 



( Vrsnak fc Skender 



20051) . the half-angle typically adding up to several degrees 



20051 ). Thus, for the CS length of several solar radii, the Petschek-CS 



thickness (the distance between the SMSs) is on the order of 100 Mm. Note that in the case 
of the vertical CS in the solar corona, where the magnetic field diverges radially, the angle 
between SMSs should increase with the height even at uniform inflow Mach number (see 
Appendix), so the CS thickness should be somewhat larger than in the plane-symmetric 
case. 

Further important characteristics of the post-CME CS that determine how it should 
look in coronagraphic images, are the horizontal length and the orientation with respect to 
the line of sight (Fig. [2]). In this respect, it is important to note that the angle 9 between 
the plane of CS and the line of sight (LoS) has to be rather small, otherwise the CS would 
not appear as a distinct feature in coronagraphic images. This causes a "selection effect" , 
which can explain the fact that post- CME rays can be recognized only in a relatively small 



fraction of events (Webb et al 



20031) . 



In Fig. 2a we depict how the effective thickness D of the ray depends on the CS thickness 
2d, horizontal length L, and the angle 0. On average, a ray should be thinner and brighter 
for smaller 9 (larger LoS column-length A). In Fig. ^Bp we illustrate ho w the variation o f 
9 along the LoS can lead to a complex appearance of the ray (see also 



Saez et al. 



2007 ). 



i.e., a wavy CS would be seen by observer as a "multi-ray" structure. Note that the angle 
9, or its variation along the LoS, can vary in time and height, so the morphology of the 
structure can change in time and can be different at different heights. 



3. Observations 



3.1. SoHO/LASCO mass images 

For the analysis of the morphology and the column-density of post-CME rays we employ 



Brueckner et al 



data provided by the Large Angle Spectroscopic Coronagraph (LASCO; 
19951) on board SoHO. We focus on the LASC0-C2 images, covering the radial distance 
range from 2 to 6 solar radii (r©). In particular, we use the so-called mass images that are 
derived from LASCO base-difference images calibrated in units of solar brightness. These 
images show changes of the column-mass along the line of sigh t under the as s umption that 
the structure lies in the plane o f sky (for details we refer to 



Billings 



1981 



Vourlidas et al 



2000 



1966 : 



Poland et al 



20021 ): bright pixels (positive values) represent areas where the 



column-mass is increased, whereas dark pixels (negative values) show depleted regions. 

The initial sample of post-CME rays considered for the analysis consisted of eleven 
events. We emphasize that these rays were spotted while analyzing the associated CMEs 
for other purposes, i.e., the rays were chosen rather randomly. Thus, no attempt was made 
to analyze the occurreiice rat e of post-CME rays (for the criteria and statistical background 
we refer to I Webb et al.ll2003r ). Nevertheless, it should be noted that ray- like features in the 



6 Vrsnak et al.: Morphology and density of post-CME current sheets 

wake of CMEs are not uncommon, particularly bearing in mind specific circumstances 
under which a post-CME current sheet could be observed. That is, the level of post-CME 
activity should not be too high, current sheet has to be oriented at small angle with respect 
to the line of sight, it should be relatively stable, etc. 

Inspecting the selected sample we noticed intermittent activity in most of the rays, 
manifested as changes of their shape, contrast, inclination, etc. These changes were either 
caused by coronal disturbances coming from remote eruptions, or were revealing internal 
activity, such as blob formation/ejection and/or outward propagating wave-like perturba- 
tions. 

For a detailed analysis we chose three post-CME rays, observed on 8-9 January 2002, 
18 November 2003, and 26 June 2005. The events of 18 November 2003 and 26 June 2005 
were selected as examples of relatively stable rays. In both events we analyzed only one ray 
image (12:50 UT and 06:30 UT, respectively). On the other hand, the event of 8-9 January 
2002 was showing a significant activity, similar to that in the remaining eight events. This 
event was selected to quantify the level of variability in post-CME rays, so we performed 
measurements at three different times (00:06, 06:06, and 12:06 UT). Note that rays in the 
selected exposures did not show significant blob-like features or similar inhomogeneities, 
although they might have been present in earlier and/or later times. 

3.2. The ray morphology 

In Fig. [3] we illustrate the variability of the ray morphology by showing two images of the 
ray that became recognizable a few hours after the east-limb CME eruption on 8 January 
2002 (first appearance in LASC0-C2 at 17:54 UT). In Fig. [3^ we first show a simple ray 
pattern that was recorded half a day after the eruption. Four hours earlier, the ray had a 
complex form consisting of three radial substructures (Fig. [HJj); this "multi-ray" pattern 
could be a result of a wavy current sheet, as illustrated in Fig.[2lD. Morphological changes of 
this kind, developing on the time scale of hours are a common characteristic of post-CME 
rays, since they were present in nine out of eleven events from our initial sample. 

The estimates of the ray width and the column-mass excess are based on the LASCO- 
C2 mass- images (Fig. [H] left), which provide for each pixel the difference Am between the 
column-mass in the actual image and the reference image. In the right-hand panels of Fig. [3] 
we show two mass-difference profiles Am{y), where Am is expressed in g/pixel, while the 
abscissa, corresponding to the y coordinate of LASCO images, is represented in pixels. The 
profiles shown were measured in the LASC0-C2 difference-images along the bin marked by 
white lines in the left panel of Fig. [31 i.e., roughly perpendicular to the ray direction. The 
bin is 10 pixels wide, and the profile shows Am averaged over the bin width. Pixels with 
Am — represent areas where the column-mass remained the same as in the reference 
image. Note that the mass astride the rays is depleted (Am < 0), probably due to the 
CME-associated coronal expansion, usually seen as coronal dimming when the eruption is 
launched from regions on the solar disc. The depletion might be partly caused also by the 
presumed inflow into the current sheet. 
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Fig. 3. a) Simple post-CME ray (9 January 2002 at 06:06 UT); b) Complex ray consisting 
of three radial substructures (9 January 2002 at 02:30 UT). The ray is marked by white 
arrow. The scale of the LASC0-C2 images is given in pixels (1 pixel corresponds to 11.9 
arcsec). The column- mass profiles measured along the bin marked by white lines are shown 
in the right-hand panels (abscissa scale in pixels). Estimates of the ray width, D, and 
column-mass excess, Sm, are indicated in a). Dashed arrows indicate the lower and upper 
limit of D. 

Given the working hypothesis, where the ray is considered to be the reconnection out- 
flow jet of a large-scale Petschek-like bifurcated current sheet, we are interested in the 
column-mass excess of the ray, Sm, defined as the difference between the peak in the ray 
column-mass Am,, and the column-mass in the adjacent corona AttIc. In other words, Sm 
determines the difference of column-mass associated with the current sheet and the inflow 
region {Sm — Amr — Amc', see the right-hand panel of Fig. [3^). The ray width D is mea- 
sured as the full width of the ray structure in the mass-difference proflle (Fig. [3^ right). The 
main source of error in estimating Sm and D is in determining the value Am^ since this 
value is frequently different on opposite sides of the ray (Fig. [3^ right). Thus, we estimated 
the upper and lower limit (dashed arrows in Fig. [3^ right), and these measurements gave 
us the error bars plotted in Fig. [H 

The measurements illustrated in Fig. [3] were performed in five images of post-CME 
rays over the LASC0-C2 range of heights to estimate radial dependence of relevant ray 
parameters. The measured rays are marked in the mass-images shown in the left column 
of Fig. [1] (hereinafter denoted as rays a - e) . In the second column of Fig. 2] the ray widths, 
expressed in units of the solar radius tq , are presented as a function of the radial distance, 
D(R), where R — r/r©. The ray in Fig. 3^ is characterized by "multi-ray" structure, so we 
measured only the central feature, which was the most prominent element of the structure. 

Inspecting the graphs one finds that in the distance range R — 2-2.5 the widths span 
between 0.1 and 0.4 vq (mean 0.26). In the range R = 5-6 the widths increase to 0.4- 
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Fig. 4. The five measured post-CME rays (indicated by arrows in the full-resolution 
LASC0-C2 images shown in the first column): a)-c) 9 January 2002 at 00:06, 06:06, and 
12:06 UT, respectively; d) 18 November 2003 at 12:50 UT; e) 26 June 2005 at 06:30 UT. 
In the second, third, and fourth column we show the ray width, column-mass excess, and 
the density excess, as a function of the radial distance. The widths and radial distances 
are expressed in units of the solar radius. The least-square fits are shown in the insets, 
together with the correlation coefficient C. 
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0.9 r (7) (mean . 7). Th e widths we found at i? ^ 2 are more than twice larger than found 



bv lWebb et al.l (|2003l ) who analyzed the coronagraph data on board the Solar Maximum 
Mission (SMM): converting their values of the width expressed in degrees, into units of 
the solar radius a t the corresponding heights, we find that the mean width reported by 



Webb et al 



(|2003l . see Tabs. 2 and 3 therein) is 0.09 in this height range. The difference 



may be related to the lower sensitivity of SMM. 
3.3. The ray density 

In the third column of Fig. [5] we show the radial dependencies of column-mass excesses 
6m(R), expressed in g/pixel. Inspecting the graphs Sm(R) we find that the slopes of the 
fitted power-law functions range from —2.1 to —3.2. The average value adds up to — 2.6±0.5. 

The column-mass excess Sm was converted to the number-density excess 6n by assuming 
that, due to the spherical geometry, the column length increases proportionally with the 
radial distance, A oc _R. In particular, we used Aq = 100 Mm at R — Rq = 2.16 (the lowest 
height of measurements) . The outcome is presented in log-log graphs in the fourth column 
of Fig. m together with the power law fits. Note that some other choice of Aq would shift 
vertically the 6n{R) dependence, without changing the slope in log-log graph. The slope 
would change only if some other radial dependence X{R) would be applied, e.g., the angle 
9 might be height-dependent. 

Inspecting the graphs 6n{R) in Fig. |4]we find that the slopes of the fitted power-law 
functions range from —3.1 to —4.2, the average value adding up —3.6 ± 0.5. [We pay 
special attention to the slope of the fit, since it does not depend on the presumed value 
of Ao, so it can be directly compared with the model results (Sect. 4.1).] In this respect 
it is important to note that in the case of the smallest slope in the sample, the ray width 
becomes approximately constant beyond R ~ 3, indicating a possible decrease of 9 with 
the height, i.e., the increase of A. This would imply that in this case the densities are 
overestimated at larger heights, since 5n cx Sm/X, i.e., the true slope of Sn{R) is probably 
somewhat steeper. 

Since the evaluation of the density excess depends on the choice of the column length 
A, in Fig. [5] we present an estimate of the number of electrons per unit length of the ray, 
which is a parameter independent of A. If we approximate the mass-difference profiles across 
the ray by triangular profiles, the mass-excess per unit length of the profile amounts to 
M = DSm/2, from which we evaluate the electron-number excess per unit length of the 
ray, N^. 

Inspecting Fig. [5^, where Ne{R) is presented separately for all five rays from Fig. [H we 
find that Ne{R) either decreases (rays [a], [b], and [d]) or stays roughly constant (rays [c] 
an d [e]). When averaged, the dat a from Fig. [5}3 show a dependence similar to that obtained 



bv lCiaravella &: Ravmondl (|2008r ). who found that after a decrease at low heights, the value 
of Ne becomes approximately constant. Comparing their Fig. 8 with our Fig. [SJs we find 
that the numbers and the trends are similar. However, the difference is the radial distance 



at which Ne{R) becomes constant: according to results bv ICiaravella fc Ravmondl (|2008| ) 



this happens at i? « 2, whereas in our case the transition occurs at i? « 4. In Fig. [S]3 we 
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Fig. 5. Number of electrons per unit length of the ray, Ne{R): a) presented separately for 
the five rays shown in Fig. [5] (measurements are labeled in the legend, following the sequence 
in Fig. |4]); b) averaged over the five rays (error bars represent standard deviations), shown 
together with the power-law fit. 

present also the power-law fit, showing that the overall trend could be described as well by 



3.4. Coronal flows 

Coronal regions beneath CMEs are generally characterized by highly dynamical intricate 
structures, revealing complex magnetoplasma flows and waves. Consequently, it is difficult 
to find persistent systematic flow patterns, especially in periods of increased solar activity, 
when disturbances from other CMEs affect the coronal region of interest. This also holds 
for flow patterns associated with post-CME rays, which show permanent morphological 
changes and are strongly affected by disturbances from distant CMEs. 

Yet, in certain situations some characteristic flows could be identified. Most frequently, 
outward-moving inhomogenei ties along the ray are observed , usually having velocities of 



Ko et al 



2003 



Lin et al 



20051 ). Such motions are often in- 



several hundred kms ^ (e.g. 
terpreted in terms of the reconnection outflow, which may or may not be characterized by 



Barta et al. 



20081 . and references therein). 



the Alfvenic speed (for a discussion we refer to 
In this respect we note that sometimes signatures of shrinking loops are also observed in 



the wake of CMEs fe.g.. lSheelev &: Wang 



2002 



Sheelev et al 



2004 : 



Sheelev fc Wang 



2007 



and references therein), consistent with decelerated reconnection downflows. 
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Fig. 6. a) LASC0-C2 running-difference images of the west-limb ray of 8 January 2002 
(13:54 UT), revealing the inflow into the ray. The shce used in the stack-plot shown in b) 
is outhned by white rectangle centered at x = -1-3.34 tq. b) Stack-plot composed of slices 
taken from successive running-difference images; x-axis represents the ordinal number of 
the LASCO image for 8 January 2002. 



Compared to reconnection outflows, detecting sig natures of reconn e ction inflows is much 
more difficult. A possible example was reported bv lYokovama et al. ( 2001), based o n the 

( 20051) . who 



Lin et al 



observations by the EIT/SoHO. Another example was presented by 
employed data from EIT and UltraVio let Coronagraph Spectrometer (UVCS) on board 
SoHO. Recentlv. lBemporad et al.l (j2008r ) reported a "side-reconncction" in the aftermath of 
a CME, induced by the CME expansion: from the detection of two converging reconnecting 
features at 1.7 tq (probably the CME flank and the streamer bound ary) the authors 
inferre d an inflow speed of 3-4 km/s, close to the 5 km/s derived bv lYokovama et al. 
(20011). 

In Fig. [S] we present a stack-plot showing an example of the reconnection-inflow pattern 
observed above the west limb on 8 January 2002 in the LASC0-C2 field-of-view. Note that 
this post-CME feature is not the ray of 8-9 January 2002 presented in Figs. [3] and IH which 
was on the opposite side of the solar disc. This west-limb event was not observed by UVCS, 
i.e., we do not know if it was characterized by the Fexviii emission, which is one of criteria 
for identifying post-CME current sheets. It should be stressed that we could not find a 
clear/measurable signature of reconnection inflow in the three post-CME rays analyzed 
in Sects. 3.2. and 3.3. In these events no suitable inhomogencities could be identified to 
provide tracing of plasma motion in the inflow region (to be discussed in Sect. 5). 

The west-limb ray presented in Fig. [SI formed after a faint ejection first observed in 
LASC0-C2 at 02:54 UT on 8 January 2002 (not hsted in the LASCO CME catalog). Each 
stack represents the cut along a line parallel to the y-axis of a given LASC0-C2 512 x 512- 
pixel running-difference image, showing the pixel-intensities in the range x — 3.30- 3.38 tq. 
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Note that the ray was oriented very close to the x-axis of LASCO images, i.e., the stacks 
are roughly perpendicular to the ray direction. The stack-plot contains all images recorded 
on 8 January 2002 (denoted by ordinal number at the a;— axis of the stack-plot; the time 
cadence is on average 24 min). 

The black/white stripes converging towards the ray location, located at ?/ ~ 400" ^ 
O.41r0, clearly show inflows into the ray. The corresponding velocities in the left part of 
the stack-plot are found in the range 15-25 kms~^ (mean 19 ±4 k ms~^) . Such velocities 
are several times higher than 5 kms~^ reported by Yokoyama et al. (|2001l) . After the slice 
No. 50 (18:54 UT), the stack-plot reveals effects of a large-scale perturbation from another 
CME, which occurred above the E-limb. Note that the "push" caused by the perturbation 
increased the inflow speed, now ranging from 25 to 30 kms^^ (mean 27 ± 5 kms^^). 

Finally, it should be noted that we also obs erved intermittent d ownflows in the form 
of shrinking loops (similar to that reported by ISheelev et al.l 120041 ) , starting from R ~ 



3. At larger heights, outflo ws could be recognized, sometimes showi n g a pattern like i n 
disconnection events (e.g. 



Webb fc Oliver 



1995 



Simnett et al 



1997 



Wang et al 



19991) . 



This implies that the diffusion region in this event was located around R 
higher than in the events shown in Fig. 3) 



3, i.e., much 



4. Comparison with the model and previous studies 

4.1. Comparison with the model results 

In Fig. [7^ we present a comparison of measurements presented in 4th column of Fig. [5] 
with the model-dependencies Sn{R) derived in the Appendix. In particular, we present the 
model results based on the isothermal (T — 10^ K) |Parkeii (|l958r ) solar wind model (see also 



Mann et al 



1999al ). For th e magnetic field i 3(i?) we take the empirical coronal magnetic 



field scaling estabhshed bv lDulk &: McLeanI (|l978l ). Note that data [a], [b], [d], and [e] are 
slightly shifted in i?-coordinate (symmetrically with respect to [c]) to avoid overlapping. 
Figure [7^ clearly shows that the observations fit well to the model curves (as stated in 
Sect. 3.3 the choice of the column length Aq only shifts the values up or down, but does 
not change the slope). Comparing the trend of the data we see that the slope corresponds 
much better to the model slopes for the CSs than for the ambient corona. Furthermore one 
finds that the CS densities are more than one order of magnitude larger than that in the 
ambient corona. 

In Fig. [7)3 we present the measured average density excess 5n{R) (circles with error bars 
represent the mean values from Fig. [7^ with the associated standard deviations) and the 
corresponding power-law fit, compared with the CS model dependencies. The observational 
data show the dependence R~^^^, whereas in the same height range the CS model data 
show ~ ° -i?"'^ -^ and the quiet-corona model behaves as R"'^-^ (see Appendix). 



4.2. Ray densities from previous studies 

Post-CME rays were studied in a number of papers, employing various techniques to esti- 
mate their density. In Table[T]we present an overview of the spectrographic results obtained 
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Fig. 7. a) Comparison of the density-excess measurements presented in the right hand 
column of Fig. [5] with the model results (thin black lines labeled by the heliocentric distance 
of the diffusion region, R^). Measurements are labeled in the legend, following the sequence 



in Fig , m The thick-gray line shows the density model of a quiet 10^ K corona (jMann et al 



1999al ). b) The mean values of the data presented in (a), shown by circles with error bars and 
the corresponding power-law fit (written in the inset), compared with the model re sults and 

( 20031): tr iangle 



t he values from pr e vious studies rev iewed in Sect. 4. 2 (diamond 



Ciaravella et al 



|2002h : squares 



Ko et al 



( 2003 ): cross ■ 



Ko et al 



Bemporad et al 



toom . 



from Ultrav iolet (UV) spectra of fiv e post-CME rays ( dates are given in the h e ader row 
report e d by 



Ciaravella et al 



(|2002h 



Ko et al 



(2003L 



Ravmond et all (|2003l ). 



Lee et al 



(|2006r ). iBemporad et al.l (|2006f ). and lCiaravella &: Ravmondl (j2008l ). In the last column we 
also present the unpublished results by Schettino et al. (2008, to be submitted). 

All rays listed in Table [1] were characterized by the presence of the narrow Fe XVIII 
emission. UV spectra of rays have bee n obtained with the UltraViolet Coronagraph 
Spectrometer (UVCS; iKohl et al.l llQQST ) on board SoHO. The spectrometer observes 



through a narrow entrance slit 42' long and up to 84" wide. The slit can be located at 
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Table 1. Spectrographic results on the UV current sheets. 





1998 Mar 23^ 


2002 Jan 8^ 


2002 Apr 21^ 


2002 Nov 26" 


2003 Nov 4^ 


2003 June 2^ 


CME Ist-LASCO 


09:33:36 


17:54:05 


01:27:20 


17:06:21 


19:54:05 


08:54:05 


CME speed (kms"^) 


403 


1794 


2393 


479 


2657 


980 


Flare 


N 


Ybl* 


X1.5 


N 


X30.8 


M3.9 


AR No. 





9782/85 


9906 





10486 


10365 


LASCO CS 


N 


2 days 


N 


N 


~ 20h 


6.5 h 


CS obs. start 


16:00:30 


17:48/JanlO 


00:45:34 


18:39:15 


20:03:50 


09:40:44 


CS obs. height 


1.51 (-1.38) 


1.53 


1.62 


1.61 


1.69 


1.68 


CS FWHM width (r©) 


0.1 


0.2 


0.3 


0.4 


0.2-0.1 


0.3 


CS full width (r©) 


0.4 


0.4 


0.6 


0.8 


0.5-0.2 


0.35 


CS LoS depth (re) 


0.06 


0.2 




0.5 


0.1 


0.1 


CS obs. duration 


20 h 


10 h 


14 min 


2.3 days 


~ 17h 


6.3 h 


CS temp. (10'' K) 


5 


3-4 


5 


8-3 


8-4 


7 


CS density(10''cm-3) 


5-10 


< 4 




6.5-7.5 


7-10 


1 



Ciaravella et al 



Ko et al.l (20031') 



(2003); 



Raymond et al 



Bemporad et al.l (120061) 



( 20031) : 



Lee et al. 



( 20061 ): 



Ciaravella fc Raymond! (|200S 



Schettino et al. (2008, to be submitted) 
* Ybl = Yes (behind the limb) 



any polar angle and at heliocentric distances from 1.5 up to IQtq. UVCS can image spec- 
tra of the solar corona in the range 945-1270 A (473-635 A in the second order). 

The first five rows of Table 1 specify the associated CME/flare events: the first appear- 
ance of the CME in the LASC0-C2, the CME mean speed, the soft X-ray importance of 
the flare if observed, the active region label, and duration of the LASCO-ray if observed, 
respectively. 

The next eight rows concern the ray characteristics as detected in the UV spectra: the 
time of the first detection of the ray by UVCS, the height at which the spectrograph slit 
intersects the ray, the ray width along the UVCS slit expressed in units of tq (both the 
full- width at half-maximum and the full width are presented), the estimate of the LoS 
depth, duration of the ray UVCS observations, and finally, in last two rows we present the 
ray temperature and density. 

The ray UVCS densities are presented in Fig. [TJd, together with the density-excess 
measureme nts presen t ed in Sects. 3.3 and 4.1. The estimate based on the UVCS spectra 



reported bv iKo et al.l (|2003r ) is indicated by the diamond symbol. The column- length was 
assumed therein to be 140 Mm which was taken to be the same as the FWHM of the spatial 
profile of the FeXVIII A974 emission across the post-CME rays. Such a column-length is 
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about 2 times larger than the value we would get for this height by the scaling that was 
applied in Sect. 3. Thus, to adjust the values to our Aq the density has to be multiplied by 
a factor of 2 (indicated by vertical bar attached to the diamond symbol in Fig. [TJ)). 



The density reported by 
is shown by the black triangle 



Ciaravella et al 



(2002), also estimated from the UVCS data. 



Ciaravella et al 



(120021 ) estimated the CS density at i? 1.5 



to 5- 10 X 10^ cm~'^ assuming the column length of 40 Mm. After adjusting to the column- 
length equivalent to that we used in Sect. 3 (required increase is around 20 %), we find the 
equivalent density around 6 x 10^ cm"'^. 

The density inferred by iBemporad et al.l (|2006l ) from the spectral observations of the 
face-on current sheet of 26 November 2002 is shown by the black cross. The attached 
vertical bar indicates that the density could be lower if a larger current sheet thickness 
would be assumed. The gray plus symbol represents the estimated density of the ambient 
coronal plasma. 

In Fig. [7)3 we also show ( black squares) th e density excess data based on the ray and 
coronal densities estimated bv iKo et al 



( 2003 ) employin g the LASC0-C2 observations and 
using th e Thomps o n scat tering function (|Billingslll966l ). In estimating the ray density at 
R = 4.4 



Ko et al 



(|2003f ) assumed the column-length of ~ 250 Mm, which is about 1.2 
times larger value than we used at the same height. Thus, adjusting their results to our 
column- length would result in 1.2 times lower density. 

Finally , in Fig. [7b the resu lts based on the polarization-brightness measurements per- 



formed by 



Poletto et al 



(|2008i ) are presented. The thick black line represents the density 

( 20081 ) by applying column-lengths con- 



Poletto et al 



excess 5n re-evaluated from Fig. 9 of 
sistent with those applied in Sect. 3. The ray data drawn by the thick line fit very closely 
to our measurements, whereas the "quiet corona" data (thin gray dashed line) are very 
close to the model densities. 

The power-law fit to all data (adjusted to the same Aq) shown in Fig. [Tb has the slope 
(without the Aq adjustment would add up to R^^ '^). This is quite close to the 
model slopes found in the radial distance ran ge 1.5-7 rp,, cha r acter ized by the power-law 
exponents 3.0-3.2 for i?x = 1.1-1.5 and the lDulk fc McLeanI (|l978r ) magnetic field. 



5. Discussion and conclusion 

The empirical characteristics of post-CME rays we summarize as follows: 

1. rays often show activity in the form of outflowing features, morphological changes, and 
changes of the inclination; sometimes reconnection inflows are observed; 

2. the width of rays increases with height, from ^0.1-0.3 r0 at i? 2 to ~ 0.4-0.8 tq 
aX Rr^Q. 

3. densities found in rays are at least several times (up to more than one order of magni- 
tude) larger than in the ambient corona in the considered height range; 

4. coronal regions surrounding rays are depleted; 

5. on average, the number of electrons per unit length of the ray, iVe, first decreases with 
the height, and then, at larger heights becomes approximately constant; 
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6. temperature of rays in the range R ~ 1.5-1.7 spans from 3 to 8 MK, i.e., it is several 
times larger than in the "normal" corona, with a tendency to decrease in time. 

The nature of the coronagraphic white-light ray activity is twofold. The formation of 
fine structure elements, mostly blob-like features, and their outward motion at speeds in 
the range 100-1000 kms~^, seems to be the most common internal activity. It could be 
interpreted in terms of the current sheet tearing, resulting in th e formation and outward 



ejection of plasmoids (e.g., iBarta et al.l 120081 : iRilev et al.l 120071 . and references therein) 



On the other hand, changes of the ray inclination, as well as morphological changes that 
could be attributed to changes of the ray geometry, seem to be caused by the large-scale 
magnetic field evolution in the wake of the CME. In some cases this type of activity is 
caused by perturbations coming from distant eruptions, most likely the large-amplitude 
waves or shocks. 

As stated in item 2, the width increases on average by about 0.4 over the distance of 
4r0, i.e., the width is roughly proportional to the radial distance. Following our interpre- 
tation of rays in terms of the reconnecting current sheet, the ray boundaries should outline 
the slow mode shocks. Given the model presented in the Appendix, in particular Fig. [51 it 
can be concluded that the slow mode shocks are oriented very close to the radial direction, 
i.e., that the angle denoted in Fig. [5] as (j>, is very small. Since the angle (f> (expressed in 
radians) equals to the inflow Alfven Mach number reduced by the factor 712/ ni that rep- 
resents the density jum p at the slow mode shoc k (the value of n^jnx depends on /3, and 
should be around 2; see 



Vrsnak fc Skendei 



20051 ). one finds that the inflow Mach number 



is low, most likely in the order of 0.01 or less. [For larger values, would be large enough 
to make the "super-radial" widening of the ray measurable.] This may be the reason why 
inflows like those shown in Fig. [5] are rarely observed. 

Note that the situation shown in Fig. [6] is different from the previously described "sta- 
tionary" pattern. In this case we see distinct elongated features inflowing from both sides 
towards the ray axis of symmetry. Thus, these features cannot be interpreted as CS bound- 
aries, i.e., signatures of the slow mode shocks, but rather they represent coronal density 
inhomogeneities aligned with the magnetic fieldlines. These features are inclined from the 
radial direction, i.e., show the "super-radial" orientation. When observed in pairs, they 
form a V-pattern roughly symmetric with respect to the ray axis. Thus, these substruc- 
tures most likely outline the magnetic field inflowing into the current sheet. In such a 
situation the angle (expressed in radians) between t he fieldline and the axis of symmetry is 



equal to the inflow Alfven Mach number Ma (e.g., Vrsnak fc Skender 2OO5I ). In the event 
shown in Fig. [6l we estimate this angle to a few degrees, so we can take Ma '^0.02-0.05. 
Bearing in mind that the inflow speed was estimated to ^ 20 kms^^, this would correspond 
to t he Alfven speed of 4 00- 1000 kms^^, which seems reasonable for this height range (see. 



e.g. 



Vrsnak et al 



2004D . 

The main physical characteristics of post-CME rays are their mass/density excess and 
increased temperature (Table I). This, as well as the morphology and flows, can be explained 
in terms of the reconnection outflow-jet, being the structural element of the vertical current 
sheet that forms in the wake of a CME. The measured values of the CS density and the 
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temperature provide also an estimate of the ambient coronal magnetic field, since the 
external magneti c pressure is rough ly equal to the gas pressure in the outflow region (see, 
e.g.. Appendix in 



Aurass et al 



20021) . Using the values ncs = 5 - 10 x lO'^ cm^^ ^^^^j j.^^ ^ 



5-8 MK, we find for the ambient magnetic field Be ^ 0.9 - 1.7 gauss. Suc h a magnetic field 



Dulk fc McLean 



is con sistent with the empirical scaling B = 0.5{R — established by 

(|l978l) . which gives B ~ 0.85-1.41 gauss for the radial distance range R = 1.5-1.7. 
Similarly, using the relationship between the temperature jump at S MSs and the externa l 
plasma-to- magnetic pressure ratio, T2/T1 



Aurass et al 



2002), 



1 + 0.4//3 (see Eq. 12 in 
one can estimate the value of jS in the inflow region. Applying T2/T1 < Tcs/Tc = 2-8, we 
flnd P > 0.4-0.06. Such values ar e roughly co nsistent with those expected for the active 



region corona in this height range (jGarv 



20011) . 



The results of the model, based on the quasi-stationary iPetschekl (|l964l) reconnection 
regime, are consistent with the observed density excesses. At the standing slow mode shocks, 
formed by the CS inflow, the plasma is heated, compressed, and accelerated/deflected to 
form the upward directed reconnection jet. In this way the dense plasma from low corona 
is transported to larger heights, causing the observed density excess. Model results show 
that the density can be increased by more than one order of magnitude, which is consistent 
with observations. However, it should be noted that the observed values depend on the 
assumed value of the LoS CS-depth Aq. Our results indicate that Aq is in the order of 100 
Mm, and that the plane of the observed rays is inclined at small angle with respect to the 
LoS. 

The parameter that does not depend on the LoS CS-depth, or the CS orientation, is 
the number of electrons per unit length of the CS, which we denoted in Sect. 3.3 as Ne- 
The measurements show that Ne(R) either decreases, or is approximately constant. On 
average, Ne{R) first decreases, to become a pproximately constant at Rc ~ 4. Given our 



Ciaravella fc Raymond! (j2008r ). it seems that 



observations, as well as the observations by 
Rc is different in different events. 

Bearing in mind the equation of continuity, there are two effects that determine the 
radial dependence of Ng. One is the radial dependence of the reconnection-outflow velocity, 
and another one is the contribution of the inflowing plasma. For example, if the outflow 
speed would be constant, the value of would be increasing monotonously with the height 
due to the cumulative supply of the plasma through the SMSs (more and more inflowing 
plasma joining to the outflow). On the other hand, if the inflow contribution is negligible 
(as it is at large heights), then a decelerated outflow would be associated with increasing 
Ne{R), whereas the accelerated outflow would lead to a decreasing Ne{R). 

Inspecting outcomes for various model inputs, we found that regardless on the model 
details, in the vicinity of the diffusion region Ne{R) increases steeply due to a strong effect 
of the plasma inflow. At larger heights, beyond ~0.5rQ above the diffusion region, the 
behavior of Ne{R) becomes dependent on the model input. A match with the observations 
(decreasing Ne{R) followed by N^. ^ const.) is found only for cases where the Alfven 
velocity and the solar wind speed increase in the given height range, and only if the diffusion 
region is below R ^ 1.4. Given that the power-law fit presented in Fig. [5}d quite well 
describes the Ne{R) dependence as oc R^^, would imply that the outfiow velocity 
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increases approximately as Uout cx R. However, note that the last three data-points in 
Fig.[5jD might indicate that N,, starts to increase gradually beyond i? ~ 5. Such a behavior 
is found in all of the considered model options, i.e., after the solar wind speed becomes 
approximately constant and the Alfven velocity starts decreasing, the value of N,, gradually 
increases with the radial distance. 

In this respect, it should be stressed that the model results generally depend on the 
applied coronal density and magnetic field model. On the other hand, the Ne{R) measure- 
ments are burdened by large errors in estimating the ray width D. So, since the comparison 
of the observed and the calculated Ne{R) might lead to ambiguous conclusions, it turns 
out that the most reliable model-parameter to be compared with observations is the slope 
of the ncs{R) dependence, which should be less steep than the slope of the ambient ndR), 
regardless on the model details. Comparing the calculated and observed slopes we find a 
very good correspondence. 

Finally, the model predicts a decrease of the density excess in time, due to the rise of 
the diffusion region height. Bearing in mind that this rise is slow, it is to be expected that 
the decrease of the density excess should be slow too. Indeed, we did not detect a decrease 
of the ray density in 12 hour interval analyzed in the ray of 9 January 2002. H owever, note 
that a decrease of the density was reported by ICiaravella fc Raymond (|2008r) in the event 
of 4 November 2003. From their Table 1 we find that the density decreased by ^ 40-50% 
in about 2 hours. According to our model, such a decrease would correspond to an increase 
of the diffusion region height for AR ^0.1, which would correspond to the rise speed in 
the order of 10 kms^^. On the other hand, given the accuracy of our density estimates 
(~ 50 %), we infer that in the event of 9 January 2002 the rise speed of the diffusion region 
was less than 2 kms^^. 



Appendix 



The overall geometry of the post-CME ray, assumed to be a signature of the bifurcated 
reconnecting current sheet, is shown in Fig. [TJ The diffusion region, where the magnetic 
field lines reconnect, is located at the radial distance R^. The oppositely directed field 
lines merge at the velocity Vin{R), bringing into the reconnection outfiow jet (the ray) the 
coronal plasma of the density nc{R). We assume that the merging velocity is faster than the 
local slow-mode speed, so the slow mode sho cks (SMSs ) form in between the two inflows. 



Petschekl ( 19641 ) reconnection model. 



bounding the reconnection outflow like in the 

Given the spherical geometry of the corona, the flow velocity vector in the inflow region 
can be represented as a superposition of the radial component (corresponding to the solar 
wind speed Vgw) and the "horizontal" component Wjn) (direction of the unit vector i9 shown 
in the lower-right corner of Fig. [5]), as depicted in Fig. [TJd. 

In the stationary state the relation V x E = must be satisfled, to provide 9B /dt = 0. In 
a spherical coordinate system (see the lower- right corner of Fig. [5]) this implies d{rE0)/dr = 
0, i.e., Rv i-gB = const. For the am bient coronal magnetic field B{R) we take the empirica l 



model by Dulk fc McLean ( 19781 ) and the simple analytical model by 



fsee also 



Vrsnak et al 



Mann et al 



(Il999bl) 



|2002l) ]. which then defines Vir,{R)- 
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We denote the angle between the radial direction and the slow mode shock (Fig. [H]) as 
<j){R). Taking into account the continuity equation across the SMS 



where v„ 



VA 



( Skender et al. 



2003 



Vrsnak fc Skendei 



20051 ). we find 



(1) 



(2) 



PSMS VA + V 

sw 

which defines the geometry of the SMSs. 

For the corona, l density Pc{R) we use the isothermal solar wind model by Parker (see 



Mann et al 



1999af ) using for the temperature Tc = 1, 1.5, and 2 MK, which provides also 



the corresponding solar wind speed VswiR)- This, together with the previously defined 
B{R), also defines the ambient coronal Alfven velocity va{R) and plasnia-to-magnetic 
pressure ratio P{R), both governing the jump relations at the SMS. Assuming that the 
guiding field is negligible, the density and temperature jump at SMS can be expressed as: 

PSMS 5(1 + (3) 



( Aurass et al 



2002 



Skender et al 



Pc 

TsMS 

—L 



2003; 



5(3 ' 
2 

5^' 



(3) 



(4) 



Vrsnak fc Skender 



20051 ). which determines the 



density and the temperature of the plasma inflowing into the reconnection outflow jet at a 
given R. 

At a given radial distance R the density of the reconnection outflow is determined by 
the flow carried from lower heights and by the flow carried across the SMS. Bearing in 
mind the geometry depicted in Fig. [H the continuity equation (mass conservation) can be 
written in the form: 



(5) 



where i^+i = {Li + Li+i)/2, p-" ^ = p™{Ri+i), and vf^^ = i;'"(i?j+i) are the current 
sheet length, the ambient density, and the ambient velocity at i? = Ri+i = Ri + AR/2, 
respectively. The term on the left-hand side represents the mass flow through the upper 
surface L^+i x d^+i, where L represents the length of the current sheet perpendicular to 
the plane of Fig. [8] The first term on the right-hand side represents the mass flow through 
the bottom surface Li x di and the second one the inflow across the SMS. The flow speeds 
v°^^ can be estimated as — VA{Ri) + Vsvj{Ri) and v°^\ = VA{Ri+i) + Vsw{Ri+i)- 



Bearing in mind the spherical geometry, the length L reads 

R 



L{R) = L(i?,) 



(6) 



where R^ is the radial distance of the diffusion region and L{R^) is the current sheet 
extension in the (^-direction at i?x, i-e., it represents the length of the current sheet X-line. 
The width of the current sheet at a given radial distance can be evaluated using 

, Ri+l 



di+i 



Ri 



Ai?< 



>i+i , 



(7) 



with dx = d(i?x) ~ 0, i.e., di — (j>iAR. The first term on the right-hand side of Eq. ([7]) is 
due to the spherical geometry (hereinafter we abbreviate diRi+i/Ri = d*^i, i.e., di^i — 
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d*_^^ + ARipi^i), and the second term is due to the inclination of the SMS with respect to 
the local radial direction (Fig. [H]). 

Equations ((T]-[7]) provide evaluation of the density distribution along the post CME- 
ray, Pout(-R), for various values of diffusion region height R^- In Fig . [HI we sho w the results 



Mann et al 



obtaine d using the coronal magnetic field model bv lDulk fc McLeanI (jl978f ) and 
( 1999b), where a s for th e coronal density we use the isothermal (10^ K) solar wind model 



bv lMann et all (|l999aD . 

Let us now consider the plasma temperature along the current sheet. In the first layer 
above the X-line the temperature equals the temperature Tsms evaluated by Eq. ^ at 
the radial distance R = Ri = Rx + AR/2. The plasma of this temperature moves outward 
and enters into the next layer through the area di x Li, undergoing adiabatic cooling due 
to volume expansion. In the second layer it meets with the plasma which inflows into the 
current sheet through the slow mode shock, where it was heated to TsMsiR2)- The same 
happens successively in each further level. 

Since thermal conductivity of coronal plasma is very high, the temperature can be 
consider as uniform across each layer, i.e., along the magnetic field lines connecting SMSs 
(magnetic field component B§). Thus, we can write for the temperature of a given layer: 

rudTd + TOsMsTsMS .o^ 

Jcs = : • (Sj 

md + TOSMS 

Here, nid and Td are the mass and temperature of the plasma contained in the volume 
defined by the lines denoted in Fig. [8] as di and d*^^^. Note that we have neglected the 
thermal conductivity along the CS axis of symmetry, since the heat fiow in this direction is 
reduced due to the magnetic field component B^), and the temperature gradients are small 
due to large vertical length-scale of CS. The mass nid is determined by: 

^d ^ Ai? , (9) 

whereas msMS ^md Tsms are the mass and temperature of the plasma contained in the 
"triangular" volume between the radial and the slow mode shock. The temperature Tsms 
can be evaluated by Eq. (|4]) at the radial distance R = Ri+i = Ri + AR/2, whereas the 
mass mgMs is determined by: 

-sl^s = Pst^sAi? . (10) 

where Pg"^s is evaluated by Eq. ^ at the radial distance R — Ri+i — Ri + AR/2. 

The estimate of the temperature Tcs in a given layer still requires evaluation of the 

temperature Td, where we have to take into account the adiabatic expansion from the layer 

"«" to "i + 1". The plasma contained in the volume defined by lines and d^: 

. ^ Li -\- Lj—-\ di -\- di—-\ , , 

V, = AR ^ 2 ' 

expands into the volume: 



Li+i -f Li d 
2 

The time At can be expressed as 



= AR + (i.+irf*+i<ri - Ud^vr)At . (12) 



At=||, (13) 
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Fig. 8. Geometry of the post-CME ray model. Dash-dotted lines depict the axis of sym- 
metry, shaded area represents the slow mode shock (SMS) and dotted lines the local radial 
direction. Local coordinate system is indicated in the lower-right corner. For details see the 
text in the Appendix. 

where zJ°"\ — {v°^\ + z;°"*)/2. The difference of the first term on the right-hand-side of 
Eq. (|12p and Vi represents the expansion in ip and i9 direction, whereas the second term 
on the right-hand-side of Eq. (fT^ represents the contribution of the expansion in the 
f-direction. Bearing in mind the characteristics of the adiabatic expansion, finally we get: 



where 7 is the ratio of specific heats, for hydrogen plasma being equal to 7 — 5/3. 

The system Eqs. ((T HT^ can be solved starting from the first layer above the diffu- 
sion region, and successively calculating the values in every new layer by using the values 
obtained for previous layer. In this way we get the dependencies p{R) and T{R) for the 
current sheet plasma. 
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